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Abstract 

Wc perform a systematic study of flavor-diagonal parity- and time-reversal- violating op- 
erators of dimension six which could arise from physics beyond the SM. We begin at the 
unknown high-energy scale where these operators originate. At this scale the operators are 
constrained by gauge invariancc which has important consequences for the form of effective 
operators at lower energies. In particular for the four-quark operators. We calculate one-loop 
QCD and, when necessary, electroweak corrections to the operators and evolve them down 
to the electroweak scale and subsequently to hadronic scales. We find that for most oper- 
ators QCD corrections are not particularly significant. We derive a set of operators at low 
energy which is expected to dominate hadronic and nuclear EDMs due to physics beyond the 
SM and obtain quantitative relations between these operators and the original dimension-six 
operators at the high-energy scale. We use the limit on the neutron EDM to set bounds on 
the dimension-six operators. 



1 Introduction 



The search for flavor-diagonal CP violation has put stringent bounds on possible CP-violating 
physics beyond the Standard Model (SM). Electric dipole moments (EDMs) of particles, nu- 
cleons, nuclei, atoms, and molecules, which break parity (P) and time-reversal invariance T 
and, by the CPT theorem, CP invariance, are the archetypical observables being looked for. 
However, despite impressive measurements, so far without success. EDMs are such good probes 
of physics beyond the SM because, at the current experimental accuracy, the known source of 
CP violation in the SM, the phase in the quark-mixing matrix, provides a negligible background 
[1]. EDMs are, however, not "SM-free" probes of new physics because of the existence of a P- 
and T- violating {f^) interaction in QCD [2j. This interaction, parametrized by an angle 9, is 
flavor diagonal and generates in principle large hadronic and nuclear EDMs |3j such that the 
null- measurement of the neutron EDM forces 9 < 10"^*^ [4j. This extreme suppression begs for 
a satisfying explanation which is currently lacking although some solutions have been proposed, 
for example the Peccei-Quinn mechanism [5j. The smallness of 9 leaves room for sources 
from physics beyond the SM. These sources are the focus of this article while the 9 term will 
not be considered. 

The search for hadronic and nuclear EDMs has grown into an active field of experiments 
on many different systems |6J. Ongoing experiments on the neutron EDM aim to improve the 
upper bound by one or two orders of magnitude [7J . Similar progress is expected for the EDMs of 
diamagnetic atoms such as ^^^Hg (which currently puts the strongest limit on the proton EDM 
[S]) and several isotopes of Rn and Ra. Furthermore, there are plans to measure the EDMs 
of the proton, deuteron, and helion (the nucleus of ^He) directly in dedicated storage rings [9j. 
These experiments have an expected accuracy which exceeds the current neutron EDM limit 
by two to three orders of magnitude. Experiments on all these different systems are needed in 
order to identify the fundamental source. A single measurement can always be reproduced 
by the 9 term or any new hadronic source. 

While EDM experiments are taking place at the low-energy frontier, other tests of the SM 
are being performed at the highest scales at the LHC. Bounds on parameters appearing in 
various SM extensions are being set and are constantly improving. In this work we wish to 
set up a framework in which these complementary tests can be compared. An ideal method 
to build such a framework is the use of effective field theory (EFT) and renormalization-group 
equations (RGEs). We assume that new physics appears at a scale considerably higher than the 
electroweak scale. Just below the former scale, after integrating out the non-SM fields, effects of 
new physics can be described by effective higher-dimensional operators consisting of SM fields 
and obeying the SM Lorentz and gauge symmetry. In this way we parametrize our ignorance 
of the high-energy theory. In principle the list of effective operators is infinite, but they can be 
ordered by their dimension. The higher the dimension of the operator, the more suppressed its 
low-energy effects are. The first operators relevant for hadronic and nuclear EDMs appear at 
dimension six {10 ^1111 [T2| ll3 | ll4j. We will not consider (semi-)leptonic dimension-six operators. 

At high energies the dimension-six operators can be measured or bounded directly, while at 
lower energies, apart from integrating out heavy SM fields, it is necessary to take into account 
the evolution of the coupling constants and possible mixing of the dimension-six operators. Here 
these effects are studied simultaneously. We begin right below the unknown scale of new physics 
and add the dimension-six operators relevant to hadronic and nuclear EDMs. Because the set 
of operators is large we focus on operators involving first-generation quarks only, leaving the 
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treatment of operators involving heavier quarks to future work. We do consider operators 
containing heavy gauge and Higgs fields which, at low energies, contribute to interactions 
among light fields. The operators are evolved, considering one-loop QCD and, when necessary, 
electroweak corrections, down to the electroweak scale where the SU{2)l symmetry is sponta- 
neously broken. At the electroweak scale the heavy SM fields decouple from theory and can 
be integrated out. We do so and match to an EFT consisting of light fields only. We evolve 
the resulting operators to a scale around one GeV where QCD becomes nonperturbative (the 
QCD scale). The QCD running to lower energies has for many operators been calculated in the 
literature |T5l [HI [TTl ISl [T9| . Here we collect, rederive, and, where necessary, supplement these 
results. We do the same for electroweak corrections |12j . 

Our approach differs from previous studies mainly in the focus on gauge symmetry and the 
systematic treatment of the various dimension-six operators relevant for hadronic and nuclear 
EDM experiments. By forcing the SM gauge symmetries onto the effective operators it becomes 
possible to elegantly derive the form of the low-energy operators. We will show, for example, 
that gauge symmetry ensures that only a subset of all possible four-quark operators among 
light quarks [18l [19] will be important at low energies. Additionally, because the set of effective 
operators considered at the high-energy scale is very general, the framework is applicable to 
many models of new physics. Once the matching to the effective operators has been performed, 
our results can be used to find the restrictions imposed by the EDM limits on the particular 
high-energy model under investigation. It is also possible to compare these low-energy bounds 
with direct bounds obtained from, for example, LHC experiments. 

In order to use the hadronic and nuclear EDM limits it is necessary to calculate these quanti- 
ties in terms of the dimension-six operators at the QCD scale. These calculations are problematic 
since they involve nonperturbative QCD but reasonable values are obtained by using QCD sum 
rules [H [20] , chiral perturbation theory [181 [ZD [SH [231 1211 [25] , or naive dimensional analysis 
[261 EZ] ■ Recently, a strategy has been proposed to disentangle the 9 term and dimension-six 
sources from measurements of the EDMs of the nucleon and light nuclei |21 1 1241 128] . This strat- 
egy depends on the particular set of dimension-six operators expected to be important at low 
energies. This set of operators is derived here (a qualitative derivation was given in Ref. [25]) 
and can be used as the starting point for calculations of low-energy hadronic observables 
such as EDMs or scattering observables [29j. 

Lattice QCD could provide more reliable calculations than chiral techniques and progress 
has been made in evaluating the nucleon EDM originating in the 9 term ^30j. Other efforts 
have focused on the contributions to the nucleon EDM arising from some of the dimension-six 
operators |31j and, hopefully, from the others in future calculations. Since lattice calculations 
are difficult and expensive we suggest to focus on the set of operators obtained here. 

The outline of the article is as follows. In Sec. [2] we list the dimension-six operators right 
below the unknown scale of new physics. The low-energy effects of some of these operators will 
be suppressed by small SM factors such as weak gauge couplings or the ratio of quark masses 
to the electroweak scale. In Sec. [3] we focus on operators without such suppression and evolve 
them to the QCD scale, integrating out heavy SM fields in the process. Once the operators are 
at the QCD scale, we use the strong limit on the neutron EDM to bound the coupling constants 
of the dimension-six operators. In Sec. |4]we repeat these steps for the operators which do suffer 
from additional suppression. The limits on these operators will therefore be weaker. We discuss, 
summarize, and conclude in Sec. [5] 
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2 Dimension-six operators 



In this article we assume that physics beyond the SM appears at a scale which, considering 
the great success of the SM, lies considerably higher than the electroweak scale Mew ~ 100 GeV. 
This assumption implies that no new particles exist with masses around or below the electroweak 
scale. Just below -Mf , the effects of new physics can be described by effective higher-dimensional 
interactions consisting of SM fields only. The interactions result from integrating out the heav- 
ier fields appearing in the fundamental theory above M^. These non-renormalizable effective 
operators have to obey the SM gauge symmetries, otherwise it it hard to understand why the 
explicit symmetry-breaking does not appear at low energies [T2j. Around MeW) SUl{2) x ?7y(l) 
is spontaneously broken down to U{1) and one is tempted to start the analysis at lower energies 
by constructing all operators that obey SUdS) x U{1) invariance. The latter strategy does not 
exploit all available information, as the breaking of SUl{2) x C/y(l) is dictated by the SM and 
does not happen in some general way. For example, in Refs. \18\ [T9] the most general form of 
SUc{3) X C/(l) invariant four-quark interactions is studied. As will be shown here, generally, 
only a subset of these interactions actually needs to be considered while others are suppressed 
due to their particular gauge-symmetry breaking properties. 

The lowest-dimensional operators important for flavor-diagonal hadronic P and T violation 
are of dimension six. The full list of gauge-invariant dimension-six operators was constructed in 
Refs. [10\ I14j. We focus onf^ interactions involving the first generation of quarks only, leaving 
the inclusion of heavier quarks (and generation-changing operators) to future work. We first 
define the fields. The left-handed light-quark fields form a doublet of SUl{2) 

while the right-handed fields uji and dji are singlets. The field (p denotes an SUl{2) doublet of 
scalar fields 99*^, J = 1,2. It is convenient to define (p^ = e^'^ip"'*, where e^"^ is the antisymmetric 
tensor in two dimensions (e^^ = +1). We will adopt the unitarity gauge, in which 



ju_ / 



(2) 



in terms of the vacuum expectation value (vev) of the Higgs field v and the Higgs boson h. 
The gauge bosons associated with the gauge groups S'[/c(3), SUl{2), and C/y(l) are denoted by, 
respectively, G^, 11^^, and B^. The covariant derivative of matter fields is 

D^^d^- igs G^e - d W'/ - ig'YB^, (3) 

where 5^, and g' are the SUc{3), SUl{2), and Uy{1) coupling constants; and and are 
SU{3) and SU{2) generators, in the representation of the field on which the derivative acts. 
The hypercharge Y for the quark and Higgs fields is given by Y^ql) = 1/6, Y{uji) = 2/3, 
Y{dR) = -1/3 and Y{ip) = 1/2. The field strengths are 

G"^, = d.G'^^-d^G^ + gsr'^G'^Gl, (4) 
W^;. ^ d,W',-d,W; + ge'^''W^w!;, (5) 
B^,u = df,B,-d,B^, (6) 
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with /"'"^ and e*-''^ denoting the SU (3) and SU (2) structure constants. 

We now turn to the dimension-six operators. Ah couphng constants introduced below are 
proportional to M^^. We start by considering various interactions between a quark, a scalar 
boson, and a gauge boson 

V 2 "v 

-^qia^'-'iriB^, + r^T • W^,)'^dR + h.c, (7) 

V 2 ''^ 

where f "''^ and rg*^ are complex coupling constants and an explicit factor v was taken out for 
later convenience. These coupling constants scale as and because they flip the chirality 

of the quark field we assume them to be proportional to the SM Yukawa couplings. If this 
assumption breaks down in particular high-energy models the coupling constants can simply be 
rescaled. Combined with the Higgs vev, the Yukawa couplings can be traded for the light-quark 
mass. It is convenient to write Eq. ([T]) in terms of the physical gauge-boson fields. Focusing on 
the terms only we find interactions containing massless gauge bosons 

A,q(C)EDM = -^[muQudu{uia'^''j5u)eF^i, + mdQddd{dicr'^''l5d)eFfj,,y 

in terms of the photon field = CyjB^ + SwW^ and its field strength F^u = d^Ay — dyA^, 
where Sy^, (c^,) denotes sinS^ (cos0^) where 9^ is the weak mixing angle. Qu,d denotes the quark 
electric charge in units of e = —gsw = —g'cw > and rriu^d the quark mass. The terms without 
scalar fields in the first and second line of Eq. ([s]) can be interpreted as, respectively, the quark 
electric and chromo-electric dipole moment (quark EDM and CEDM) with coupling constants 

, c^ImT^ + s^Imri;^ c^Im - s^Im ~ Imf" 
du = — 7^ dd= — ^, dg = . 9 

After electroweak symmetry breaking Eq. ([T]) also contributes to interactions involving heavy 
gauge-boson fields. Interactions involving the neutral boson = CwW^ — s^Bf^, in particular 
its field strength Z^^u = d^Zy — dyZ^, are given by 

>C6,qZEDM = -| [2:umu(mCT^'^75M) + 2;dmd((ii(T''''75d) Z^,;,] ^1 + , (10) 

with coupling constants 

c^„Im^^ - s^ImT^ c^Im + s^Im 

Zu = Zd = -. (11) 

grriu gmd 
Finally, there are terms involving the charged gauge bosons Wj^ = (Wifj, =F iW2^)/V^ and their 



4 



field strengths = (W^^ =F ^'^fiv)!^^ which are of the form 



L 



6,qWEDM 
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V2 
_9_ 
\/2 



UluWu 



rridWd 



{uLia^^'dR) W+ + i^idiaf^'^-f^d) W+W' + h.c. 



V2 



1 + 

1 + 



where 



imriv 



imr^ 



(12) 



(13) 



The quark EDMs and weak EDMs are related by gauge symmetry as can be seen from Eqs. (|9|, 
and ((13]). 

The next dimension-six operator is the Weinberg operator [11] which contains only gluons 

^ dw 



-GGG 
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rabc fj.ual3 y^a f~ib /~ic 



(14) 



where ei^'^'^l^ is the totally antisymmetric symbol in four dimensions (e^^^s _ This operator 
can be interpreted as the CEDM of the gluon |32| . 
Next we consider interactions among quarks |13j 



mqq = ^iiQLUR)£ij{qLdR) + T,s{q\taUR)eij{qltadR) + h.c. 



(15) 



The terms can be rewritten as 
.Im Si 



6,qgqg 



+i 



2 

Im Sg 



{uu dj^d + uj^u dd — du uj^d — dj^u ud) 
[ut^u dj^f^d + uj^t^-u df'd - dt^u u^y^t^-d - dj^t^-u ut"-d) + (16) 



where the dots stand for PT-even terms. These are the only dimension-six four-quark 
interactions among the first generation of quarks that are allowed by the SM gauge symmetries. 
These four-quark operators are not suppressed by light-quark masses. We denote the operators 
in Eq. (16) by FQPS because of their pseudoscalar-scalar form. Based on their color-structure 
we name the first and second term in Eq. (16), respectively, the color-singlet and color-octet 
FQPS operator even though, of course, both operators are scalars under SU{3)c- 

At lower energies additional four-quark interactions will appear when the Higgs and heavy 
gauge bosons that appear in various dimension-six operators are integrated out. An example of 
such an operator is p3] 



which after electroweak-symmetry breaking becomes 



2V2 



El UR-yi^dR + h.c. 



1 + '^ 

V 



(17) 



(18) 



Below the electroweak scale, the W boson is integrated out and the terms proportional to ImHi 
will contribute to four-quark interactions with different structure than Eq. (16) |33]. The 
associated factor Gp is compensated by the factor v"^ in Eq. (18). 
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At the level of dimension-six operators additional interactions among quarks and Higgs bosons 
appear 



(19) 



The terms proportional to {v + h) can be absorbed into the Yukawa couplings appearing in the 
SM. The terms that remain are 



h{v + h) (v+'^] {ilmY"' u-f\ + ilmY"^ d-/^d] . 



(20) 



Finally, there are operators consisting purely of gauge and scalar bosons. Analogous to Eq. 

(21) 



(14) one can write an interaction among SU{2)l gauge bosons 

d 



www 



g - - ■ ■ a/J" ^ip" u ! 

which can be written in terms of physical gauge bosons as |12j 

- idu,e^"'''^W+PW-^ {s^F^, + c^Z^, - 2igW+W-) . 
Other operators involve the Higgs field 



(22) 



_/x!/a/9 //!/ 



(23) 



Without any Higgs field, i.e. using ip^ip/v^ = 1/2, the first three operators in Eq. (23) can be 
absorbed into the 5C/c(3), 5J7l(2), and (1) topological theta terms appearing in the SM. This 
is not the case for terms which contain at least one explicit Higgs boson. These are given by 



-2s^c^{g'H'B - g'e'w)F^uZo,p + 2g'e'w W^^W^ 



.2/)/ 



2/)/ 



/!(?;+-) + .., 



(24) 



where the dots stand for terms involving three gauge bosons which will not be used further on. 
The fourth term in Eq. (23) obtains nontopological terms both with and without an explicit 
Higgs boson [12] 



99 GwB 



{SwF^u + Cu]Zfj,^){CwFai3 - SwZa(5)h{v + 



-ig{W^W^ ){c^Fap - s^Zo,^){v + /i)^ 



(25) 



Linear combinations of the operators without Higgs bosons in Eqs. (22) and (25) can be inter- 
preted as the electric dipole and magnetic quadrupole moment of the boson [M] . 



The goal is now to evolve the operators in Eqs. ((TO]), (|T2]), ([14]), ([16]), ([18]), ([20]), ([24]), 
and (25) to the low-energy scales Mew and subsequently to A/qcd- In the process we have to 
integrate out the heavy fields and include the effects of QCD and, in some cases, electroweak 
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renormalization-group running. Although all above operators are proportional to two inverse 
powers of not all operators will be equally important at low energies. Due to the need to 
integrate out the heavy gauge and Higgs bosons, certain operators induce low-energy interactions 
which are, apart from the suppression, additionally suppressed by small SM quantities such 
as electroweak gauge couplings or the ratio of light-quark masses to the electroweak scale. It is 
therefore convenient to divide the operators into two sets. The first set contains the operators 
that induce low-energy interactions without suppression. After studying this set in Sec. [3j we 
will turn in Sec. |4] to the remaining operators which do suffer from additional suppression. 



3 Unsuppressed operators 

We begin at the scale M'p where all SM fields are still present in the effective field theory. At 
this scale, the operators in the first set are given by the quark EDMs and CEDMs in Eq. ([8|), the 



gluon CEDM in Eq. (14), the four-quark operators in Eq. (16), the right-handed weak current 



in Eq. (18), and finally the quark-Higgs and gluon- Higgs interactions in Eqs. (20) and (24). The 



first goal is to evolve these operators down to the electroweak scale Mew ~ where My/ is 
the mass of the W boson. At this scale we integrate out the Higgs and heavy gauge bosons. In 
principle it would be better to integrate out each field at its particular threshold, but this would 
complicate matters without gaining much. Somewhat above Mew we pass the threshold of the 
top quark which has important consequences for certain operators. 

3.1 From My to Mew 

To calculate the QCD running of the operators it is convenient to write the Lagrangian as 

U = Y.'^^d^^O.ifi), (26) 

i 

with /i the renormalization scale and Oi{fi) the set of operators 



Ow 
Ops, 

OpSs 

Owr 

OqH 
OgH 



Oq = --emgQgqa'"'-f5q F^u, 



gsTHq qa^^-i^taq 



„ysjabc^ ^af)^ up^v i 



{uu dy^d + wy^u dd — du wy^d — dy^u ud), 

[m^u dy^ed + uy^f-u (ifd - dt"-u uy^f'd - dy^f-u ut^d) , 



2^/2 



{iuR-i^'dRW'[ - idny^'uRW^ ), 



V hqiy q, 
gy^^^G^C^^^vh, 



(27) 
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where q G {u, d}. By comparison with Sec. [2] the couphng constants at are found to be 



dw{M^) 
CgH{M^) = e'{M^), 



ImSi,8(Mr) 



CqH{M^)=l^Y"i{Mr^), 



(28) 



which are all proportional to M^^. We use the following operator basis 

C{ij) = {C,{fi),Cgifi),Cw{f^),CpsM,Cpss{fi),CwAf^),CgH{fi),CgH{f^)f. (29) 
The RGE can be expressed as 



d In/i 



7C(/x), 



(30) 



where 7 is the anomalous dimension matrix which is expanded in powers of (7^/(47r)^ = Qs/(47r). 
Up to first order in as, 7 can be written as 



7 



as 
Air 



/ 7dipolc 


7mix 








7mix gH^ 





IPS 

















^Wr 

















IqH 





V 











7gH 1 



(31) 



Most of the submatrices appearing in Eq. (31 ) have been calculated in the literature. The running 



and mixing of the quark and gluon dipole operators were calculated in Refs. [HI [151 HSl fT7] 



7dipole 




-8C2(iV) 
mC2{N) - AN 





(32) 



where C2{N) = {N'^ — 1)/2N and /3o = ^(11^ — 2nj) in terms of N and nj the number of, 
respectively, colors and flavors. Here we confirm these results. 

The submatrix 'jps describes the running and mixing of the four-quark operators. In 
Refs. (18^ ti9j these matrices were calculated for a different basis of operators consisting of the 
most general flavor-diagonal four-quark operators containing up, down, and strange quarks. 
In the limiting case of only two flavors, the set contains ten different operators whereas our set, 
based on gauge invariance, consists of two interactions only. As we will see, below Mew two 
additional four-quark operators will be induced by Owri while the remaining six operators are 
in general suppressed. The suppressed terms are discussed in Sec. |4j In the basis of Eq. ( 29 ) 
the anomalous dimension matrix of the FQPS operators is 



/ (3jV+4)(jV^-l) 
\ ^ N 



/3o 



{N+l)^{N-2){N-l) 



/3o 



(33) 
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Figure 1: One-loop diagrams contributing to the quark electric and chromo-electric dipole mo- 
ments. Solid, wavy, and curly lines represent the propagation of quarks, photons, and gluons 
respectively. The square denotes a four-quark interaction, other vertices representing T 
interactions from the Standard Model. 



We have checked that Eq. (33) agrees with results obtained in Refs. |181 119j after a suitable 
basis transformation. 

It is necessary to consider mixing of the four-quark interactions with the dipole operators. 
This is described by the submatrix 7mix- The diagrams responsible for this mixing are shown in 
Fig. [T] The FQPS operators are invariant under chiral symmetry {i.e. under global SU{2)l x 
SU{2)^ rotations) j24| such that they can only contribute to (C)EDMs, which break chiral 
symmetry, through a mass insertion. For this reason, the qCEDMs do not induce the four-quark 
operators and the mixing is one way only. Up to 0{as), the FQPS operators and the gluon 
CEDM do not mix although they are both chiral invariant. With these considerations we find 



Tmix — 2 

rua 



-1 f-C2(iV)h (34) 
\ 



where q' denotes the quark flavor which is different from the flavor of the induced quark 
(C)EDM. For example, for the mixing between the up quark EDM and Opsi the relevant 
entry is (7mix)ii = 2(mrf/mu)(— 1/3)/(2/3) = — md/m„. Again, after a basis transformation. 



Eq. (34) is in agreement with results in Ref. |19j . 

The next operator is the right-handed weak current. This operator will obviously not mix 
with the other operators via QCD corrections due to the presence of the I^-boson field. At 
the level of one-loop electroweak diagrams it will induce corrections to the quark (C)EDM, but 
these corrections are smaller than tree-level diagrams contributing to four-quark operators |33j . 
These tree- level diagrams will be considered below. At 0{as) the operator Owr does not run 
and its anomalous dimension is zero, 

7Wr = 0. (35) 

Then there are the Higgs-quark interactions with anomalous dimension 

7te = -6C2{N). (36) 

At tree- and one-loop level these interactions induce four-quark operators and the q(C)EDM, 
but such contributions are suppressed by a factor niq/v. At the two-loop level it is possible 
to circumvent this suppression via a top-quark loop shown in Fig. [3] which contributes to the 
qCEDM. This diagram will be discussed below. 
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Figure 2: One- loop diagram contributing to a quark chromo-electric dipole moment. The dashed 
hne denotes the propagation of a Higgs boson. The square denotes the gluon-Higgs interac- 
tion from Eq. (27). Other notation is as in Fig. [T] For simplicity only one possible ordering is 
shown here. 




Figure 3: Two- loop diagram contributing to the quark chromo-electric dipole moment. The 
double solid line denotes a top-quark propagator. The square denotes the quark-Higgs 
interaction from Eq. (27). Other notation is as in Figs. [l]and[2j For simplicity only one possible 
ordering is shown here. 



Finally, we consider the gluon-Higgs interaction. In order to calculate its anomalous dimension 
and its mixing with the dipoles it is useful to rewrite 



which makes it easier to obtain 



l9H = Q, 7mix3/f=(0 16 0)^, (38) 

where the mixing with the quark CEDM is due to the loop diagram in Fig. [2| The vanishing of 
^gH at the one- loop level is in agreement with Refs. [351136] , 

Having calculated the anomalous dimensions of the operators, it is possible to solve the RGE 
and run the operators down to the electroweak scale Mew- During this process we pass the 
threshold of the top quark. The top quark is integrated out and we match the EFT with six 
flavors to the EFT with five flavors. As none of the operators contain the top quark the operator 
basis remains unchanged as we move from one EFT to the next. The matching conditions are 

C\mt) = CHm^), (39) 

where are the coupling constants belonging to the EFT incorporating the top quark and 
to the EFT without it and mf (m^) denotes the energy scale just above (below) the mass 
of the top quark. The running is somewhat different for the various EFTs as the anomalous 
dimension matrix depends on explicitly as well as implicitly through as{n). We employ the 
one-loop running of as for the sake of consistency. The above treatment works for all operators 
except for the Higgs-quark interactions for which, as mentioned above, the passing of the 



10 



Table 1: The size of the quark electrie and quark and gluon chromo-electric dipole moments at My/ in 
terms of various dimension-six operators at 1 and 10 TeV. A "— " indicates that there is no contribution 
at this order. 





C,(lTeV) 


(7,(1 TeV) 


(7vK(lTeV) 


(7,H(lTeV) 


CgH{neV) 


C, {Mw) 


0.80 


0.18 


-0.010 


-1.4 X lO"'' — 

TTlg 


-0.032 


{Mw) 




0.82 


-0.090 


-2.1 - 10-*^ 


-0.30 


Cw {Mw) 






0.64 








Cq(lOTeV) 


(5,(10 TeV) 


(7w(10TeV) 


C7qH(10TeV) 


C7gij(10TeV) 


C, {Mw) 


0.69 


0.26 


-0.025 


-1.6 X 10"'' — 


-0.085 


Q {Mw) 




0.72 


-0.12 


-2.4- 10"*^ 


-0.48 


Cw {Mw) 






0.46 







top-quark threshold does have important consequences. The two-loop diagram in Fig. [3] gives a 
large contribution to the qCEDM 

= -^;^/(^) (40) 

where f{-^) — 1 ^37] is a function of the top-quark and Higgs mass. Here the light quark mass 
should be evaluated at the scale mt, however, we prefer to present results in terms of quark 
masses evaluated at the scale Mqcd • The one- loop running of the quark masses gives 

mq{mt) ~ 0.55m<y(MQCD)- (41) 

Similar diagrams as Fig. [3] but with gluons replaced by electroweak gauge bosons contribute to 
the quark EDM, but these contributions are smaller by approximately an order of magnitude 
due to the smallness of a^(mj) compared to as(mj) and we neglect them. Because the two- 
loop diagram gives by far the most important contribution to low-energy P and T violation 



coming from Eq. (20), below rrit the quark-Higgs interactions can be neglected. Their effects are 



saturated by the qCEDM contribution in Eq. (40). 



We can now express the dimension-six coupling constants at the electroweak scale (Mew = 
M\y) as function of those at the scale of new physics (Mf ). Here we consider numerical solutions 
for two specific values of iW^ (Mj- = 1 TeV and = 10 TeV). The reason being that the 
analytical solution of the RGE, while easy to obtain once the anomalous dimensions are known, 
is too lengthy to write down. The following masses are used as input [38] 

Mz = 91.2 GeV, Mw = 80.4 GeV, mt{mt) = 160 GeV, (42) 



and for the value of as 1381 



a'^^'^iMz) = 0.118. (43) 



The one-loop running then gives, 

a7^^{l TeV) = 0.089, a"^"''(10 TeV) = 0.073. (44) 

We will suppress the Uf dependence of from here on since the number of flavors is always 
clear from the context. The results are shown in Tables [T] and [2| The operator Owr is not 
shown because it does not run and gives the trivial result Cvi/fl(-/V^r) = Cwui^w)- The results 
are not very dependent on the choice of M^, because runs slowly at high energies. 
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Table 2: The size of the couphng constants of various dimensions-six operators at Mw in terms of the 
four-quark operators Ops^ g at 1 and 10 TeV. A "— " indicates that there is no contribution at this order. 
The notation q' appearing in the first two rows denotes the quark flavor which is different from the flavor 
of the induced quark (C)EDM. 



Cq {Mw) 
Cw {Mw) 
Cps, {Mw) 
Cpss {Mw) 



Cpsi(lTeV) 



0.0043 - 0.034 
0.042^ 

1.2 

-0.26 



CpgJlTeV) 
0.00053 + 0.043^ 
-0.0045^ 

nig 

0.046 
0.73 



C, {Mw) 
Cq {Mw) 
Cw {Mw) 
Cps, {Mw) 
Cpss {Mw) 



CpSi(lOTeV ) 
1^0.012 -0.05l| 
0.074^ 

nig 

1.3 

-0.44 



CpSs(lOTeV) 
0.0012-^0.063^ 
-0.0052^ " 

0.077 
0.58 



3.2 From Mew to Mqcd 

The next step is to cross the electroweak scale at which we integrate out the Higgs and heavy 
gauge bosons. Because the Higgs-quark interactions have effectively disappeared from the 
operator basis around the top-quark threshold, the only operators involving heavy fields are the 
right-handed weak current and the gluon-Higgs operators. The latter contribute to the quark 
(C)EDM not only via running from to Mew but also via finite corrections at the Higgs 
threshold. These corrections, however, appear one order higher in ag and can be neglected. The 



W boson in the right-handed weak current in Eq. (18) can be efficiently integrated out via a 
tree- level diagram among light quarks ^33j. Since the SM couples the W boson to left-handed 
quarks only, the appearing effective operator is of the form 

^6,FQLR = igsClRi {uR-^^dR dLlf^UL - dR'J^UR UL^^di) + ... 

= ^-Y^LR^ {u-f^-f^d d-f^u - d-i^'r'u u-/f,d) + . . . , (45) 

where the dots stand for higher-dimensional terms appearing from expanding the 1^-boson 
propagator. The factor g'^ is introduced for later convenience. The coupling constant right 
below the electroweak scale is given by 

where Vud — 0.97 is the up-down CKM element and M^y {M^) denotes the scale just above 
(below) the VF-boson mass. As promised, the factor v'^ in Eq. (18) has cancelled against the 



M^r appearing in the propagator of the W boson via the SM relation Myy = gv/2. 

After integrating out the heavy fields the operator basis has changed. The right-handed weak 
current, the Higgs-quark, and the Higgs-gluon interactions are removed and two new four-quark 
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operators are added 



OlRs = i 9s {uR^tadR dll^taUL - dRl'^taUR ULlfitadl)- 



(47) 



The coupling constant of OlRi is given in Eq. (46) while CLRf^{Mw) = 0. The first operator in 



Eq. (47) is produced by integrating out the W boson as discussed above, while the second arises 



when OlRi is evolved to lower energies. Because these four-quark interactions couple left-handed 
to right-handed quarks they are denoted by FQLR. The new operator basis becomes 



C{f^) = {Cq{fl),Cq{n),Cw{f^),Cps^{fi),CpSf,{fJ-),CLRi{fl),CLRs{fJ')y 



(48) 



Apart from the appearance of the FQLR operators, the basis has not changed by moving through 
the M\Y threshold, thus for the remaining couplings we have, Cj(M^) = Ci{M^), i / LRi s- 
The new RGE is written as as 



dln/i 

with the new anomalous dimension matrix 



7C(m), 



(49) 




(50) 



where, apart from "yiR, all entries can be found in the previous section. The diagrams in Fig. [T] 
vanish for the FQLR operators such that they do not induce the quark (C)EDM. Furthermore, 
the FQPS and FQLR operators do not mix because of their different chiral properties [25], 



justifying the form of Eq. (50). The submatrix ^lr describes the running and mixing of OlRi 
and Olrs- It is given by 



7LR 



-/3o 



3C2{N) > 

2 ^ 

Po, 



(51) 



which is again in agreement with Refs. |18|. IT9] after a suitable basis transformation. 



We can now solve Eq. (49) and run the operators down to the energy scale Mqcd- During 
this process we pass the thresholds of the bottom and charm quarks which are handled the same 
way as the top-quark threshold. None of the -pjT operators involve these quarks explicitly such 
that the operator basis does not change between thresholds. We use [38] 

nibinib) = 4.18 GeV, mdmc) = 1.28 GeV, (52) 

and run the operators down to Mqcd = 1 GeV. At this scale the strong coupling constant is 



as(MQCD) = 0.358. 



(53) 



The results are collected in Tables [3] and |4j By combining Tables [T]|4] it is possible to read off the 
size the coupling constants at Mqcd in terms of those at M^. Although the bases in Eqs. (29) 
and (48) are useful to calculate the anomalous dimensions, the explicit appearance of factors of 
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Table 3: The size of the couphng constants of various dimension-six operators at Mqcd in terms of 
those of at Mw A "— " indicates that there is no contribution at this order. 





Cq {Mw) 


Cq {Mw) 


Cw {Mw) 


Cps, {Mw) 


Cpss (Mm/) 


Cq (Mqcd) 


0.48 


0.37 


-0.060 


(0.040 - 0.071 1;^)^ 


- (0.0019 + 0.078^) ^ 


Cq (Mqcd) 




0.53 


-0.15 


0.14^ 

Trig 


-0.0017^ 


Cw (Mqcd) 






0.26 






Cpsi (Mqcd) 








1.5 


0.13 


Cpsa (Mqcd) 








-0.71 


0.28 



Table 4: The size of the couphng constants of the FQLR operators at Mqcd in terms of those at Mw ■ 





ClrAMw) Clrs{Mw) 


ClrAMqcd) 

CLRsiMQCD) 


0.37 0.10 
0.47 0.92 



gs in the operators is unconventional. Here we give the results for the operators without these 
factors. Defining the quark (C)EDM as in Eq. ([s]) and using the boundary conditions from Eq. 
(|28]), we obtain for M^r = 1 TeV 

dq{MQCD) = 0.39 tfg(l TeV) + 0.37 dg(l TeV) -0.13 (9' (1 TeV) 
-0.072 dvK(l TeV) - (8.4 • 10"^)— Imy'«(l TeV) 

+ {0.20, 0.097} Im El (1 TeV) + {0.073, 0.069} Im ^8(1 TeV), 
dqiMqcD) = 0.88 (ig(l TeV) -0.340' (1 TeV) 

-0.29 dvF(l TeV) - (2.4 • 10"^)— Imy"'(l TeV) 

TUg 

+ {0.74, 0.17} Im Si (1 TeV) + {0.011, 0.0025} Im S8(l TeV), (54) 

where the first (second) term in brackets in front of Im Si.g corresponds to the up (down) quark. 
For M'^ = 10 TeV we find slightly different values 

^^(Mqcd) = 0.33 (ig(10 TeV) +0.41 (ig(10 TeV) - 0.22 6''(10 TeV) 

-0.089 dvF(10 TeV) - (9.5 • 10"^)— Imy"'(10TeV) 

rriq 

+ {0.30, 0.14} ImSi(10 TeV) + {0.092, 0.085} Im Sg (10 TeV), 
dqiMqcD) = 0.85 (ig(10 TeV) - 0.54 6l'(10 TeV) 

-0.30 dvK (10 TeV) - (2.7 • 10"^)— Imy"?(10TeV) 

rriq 

+ {1.1, 0.25} Im Si (10 TeV) + {0.034, 0.0078} Im S8(10 TeV). (55) 

Without looking at specific models of new physics it is difficult to compare the various contribu- 
tions to the low-energy quark (C)EDMs. Assuming the various coupling constants to be of equal 
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size at and using v/niq ~ O(IO^) one would conclude that the low-energy quark (C)EDMs 
are dominated by the Higgs-quark interactions. However, in certain models of new physics 
the Yukawa-like couplings Y'^''^ scale with the quark mass and are thus 0{mq/vM^) instead of 
0{1/Mj,) as assumed in Sec. j2jp9]. The resulting suppression makes the contributions propor- 
tional to y'"'"^ negligible. In such a scenario, the low-energy quark EDM gets contributions of 
roughly equal size from all other dimension-six operators at although the dependence on the 
gluon CEDM and color-octet FQPS is somewhat smaller. The low-energy quark CEDMs are 
mainly determined by the high-energy CEDMs and, depending on the quark flavor, the color- 
singlet FQPS. The contributions from the gluon CEDM and the gluon-Higgs interaction are also 
significant, while the influence of the color-octet FQPS can be neglected. We stress again that 
these statements are very model dependent since the coupling constants at the high-energy scale 
might, depending on the model of physics beyond the SM, possess a large hierarchy. Some may 
even be zero. In any case, it is unlikely that all coupling constants at are of similar size. 



The next operator is the gluon CEDM, which is defined without Qs as in Eq. (14). At this 
order, the gluon CEDM does not depend on other sources such that 

diy(MQCD) = 0.33dv^(lTeV), diy (Mqcd) = 0.27 diy (10 TeV) . (56) 

The gluon CEDM is somewhat suppressed at low energies in agreement with Refs. |12 | I16 | [T7]. 
Additional contributions to the gluon CEDM do appear if the CEDMs of heavy quarks are taken 
into account [TB] . 

The FQPS operators do not mix with other four-quark operators such that the solution of 
the RGE is very simple. For = 1 TeV 

Im Si (Mqcd) = 7.2ImSi(lTeV)-h0.66ImS8(lTeV), 

ImS8(MQCD) = -3.7ImSi(lTeV) +0.69ImS8(lTeV), (57) 

and for = 10 TeV 

Im Si (Mqcd) = 10 Im Si (10 TeV) 0.95 Im Eg (10 TeV), 

ImS8(MQCD) = -5.3ImSi(10TeV)-F0.53ImS8(10TeV). (58) 

Again assuming that the coupling constants are of the same order at the high-energy scale, 
the FQPS operators at Mqcd are dominated by ImSi(M^). Barring unexpected fine-tuning, 
ImSi (Mqcd) and ImS8 (Mqcd) are of approximately the same size and both operators should 
be taken into account at low energies. It is interesting that the coefficient of the color-singlet 
operator grows significantly when evolved to lower energies. 

Finally, we look at the FQLR operators. Because the operator Own^ from which the FQLR 
operators originate, does not run between M^ and Mew the results do not depend on the 
particular value of M^. (Of course, this is true only as far as the running is concerned. The 
size of the couplings is proportional to M^^). The results are given for the operators without 
explicit factors of Qs- The redefined couplings are 

ImHL/j^(/x) = gsiufCLRiif^), ImHLRg(^) = gsififCiRsil^), (59) 

for which we find 

IuiElr, (Mqcd) = -l.lKdImHi (M^) , ImEiR, (Mqcd) = -1.4KdImHi (M;r) • (60) 
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The conclusion is that both four-quark operators are of approximately the same magnitude 
at low energies even though Olrs does not get a direct contribution at the electroweak scale. 
Clearly, the coupling constants of the two FQLR operators are not independent. Both depend 
on the same fundamental parameter Hi (My). 



The right-handed weak current in Eq. (17) not only gives rise to the FQLR operator but also 
to semi-leptonic interactions by coupling it to the left-handed lepton current. The operator 
produced in this way does not contribute to hadronic EDMs but, instead, contributes to the ^ 
triple correlation ^ D J ■ {pe x p^) in nuclear /3-decay [33j. Both the FQLR operator and its 
semi-leptonic cousin arise from the same dimension-six operator and depend on the parameter 
ImSi(My). In Ref. [33j it is argued that the best limit on ImHi(My) comes from the limit on 
the neutron EDM. Furthermore, the authors of Ref. |33j used the neutron EDM limit to set a 
strong constraint on the coefficient D. However, QCD corrections were neglected and one might 
wonder whether their effect is significant. The semi-leptonic operator does not run under QCD 
RGE which makes it sufficient to consider the running of the FQLR operator. The running of 
the color-singlet FQLR gives an enhancement of 10% which does not alter the results of Ref. |33j 
in a significant way. However, the fact that at low energies not only the color-singlet but also the 
color-octet operator is present, with a larger coefficient, may be more important. It is difficult 
to say how this affects the bound obtained in Ref. [33] since this depends on nonperturbative 
physics linking the FQLR operators to the neutron EDM. In the event of a partial cancellation 
between the neutron EDM contributions from the color-singlet and color-octet FQLR operators, 
the bound on the D coefficient would be weakened. 

To summarize this section we write the Lagrangian at the scale Mqcd as 

>C6,set 1 = - {duemuQu uia'^'j^u F^y + daemaQd dio-^^^bd F^y 



+^r'"e>^'"'^GlpG%Glf 
Iin 

+i — ^ — {uu dj^d + M75U dd — du uj^d — d'y^u ud) 

{uf'u dj^t^d + u-f^f^u dfd - df'u u-i^f'd - dj^f^u ut^d) 
+i Im -.LRi {uRl^dR dil^UL - dR'y^UR ULlfidi) 

+i Im "LRs [uR-i^t^-dR dLlfjit^-UL - dR^iH^UR ULlf.f'dL) , (61) 
where the coupling constants are all evaluated at Mqcd- Their values in terms of the coupling 



constants at the high-energy scale can be read from Eqs. (54)-(60). We expect the above 
Lagrangian to capture the most important contributions from physics beyond the SM to first 
generation observables. For example, the EDMs of the nucleon and light nuclei should 



be calculated as a function of the various coupling constants appearing in Eq. (61). At lower 



energies, Eq. (61) induces/^ interactions among pions, nucleons, and heavier baryons. These 
interactions have been constructed in Ref. [25j by use of chiral perturbation theory. 

An important result of this section is the form of the four-quark operators. Although 
they look complicated, their form is very intuitive if one insists on gauge invariance at the scale 
where the dimension-six operators are generated. This is also indicated by the simple block- 
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diagonal form (apart from mixing with the quark dipoles) of the anomalous dimension matrix. 
We conclude that around a scale Mqcd ~ 1 GeV there are four four-quark operators that 
need to be taken into account with three independent couplings. We will demonstrate in Sec. 



4.1 that the remaining six combinations also appear but are, in general, suppressed. 



3.3 Bounds from the neutron EDM 

In this section we use the stringent limit on the neutron EDM, dn < 2.9 • lO^^^efm ^ to set 
bounds on the various dimension-six coupling constants. This requires the calculation of the 



neutron EDM in terms of the couplings appearing in Eq. (61) which, in turn, can be related to 
the couplings at the high-energy scale using the results obtained in the previous sections. 
Calculating the neutron EDM in terms of the quark-gluon operators is a problematic task due to 
the nonperturbative nature of QCD at low energies. Despite this difficulty, several approaches 
exist to tackle this problem (for a review, see Ref. [1]). Here we use recent results obtained in 
chiral perturbation theory (xPT) |2H 122} |2^ I25j . which has the advantage that all operators 



appearing in Eq. (61) have been treated within the same framework. Furthermore, it allows 
for the reliable calculation of light-nuclear EDMs |23t [2H [28] which have become the subject of 
experimental investigation [9j. 

In xPT, the low-energy EFT of QCD, the effective degrees of freedom are pions and nucleons 
(and heavier baryons) whose interactions are determined by the symmetries of QCD and how 
they are (spontaneously and explicitly) broken. Pions are interpreted as the Goldstone bosons 
of the spontaneously broken chiral symmetry of QCD, 5C/i(2) x SUr{2). For a review of 
see, for example, Refs. |40[l41j. The extension of xPT to include the effects of the dimension-six 
operators has been performed in Ref. |25) (the 9 term has been studied in Ref. [l^ ) . At leading 
order in the calculation of the neutron EDM only three hadronic interactions play a role 
|21l [25] for all dimension-six sources in Eq. ( |61[ ) . These are given by 

£^ = - ^Nt ■7vN-2N{do + diTs) S^N v^F^^, (62) 
in terms of the nucleon doublet N = (p n)^, the pion triplet vr, and the pion decay constant 



-Ftt = 186 MeV. In Eq. (62) the heavy-baryon framework [431 141] has been applied where, instead 



of gamma matrices, it is the nucleon velocity v'^ and spin that appear. The first interaction in 



Eq. (62) is a chiral-symmetry-breaking vrA^ interaction, while the other two are short-range, 
i.e. due to dynamics of shorter range than pions, contributions to the isoscalar and isovector 
nucleon EDM respectively. Which of the hadronic interactions plays a role in the calculation of 
the neutron EDM depends on the particular dimension-six operator under investigation. 

The actual leading-order calculation of the neutron EDM is fairly straightforward and gives 

dn = ido-d,) + j^,ln-^, (63) 

where gA — 1-27 is the strong pion-nucleon coupling constant and m.,^ ~ 137 MeV (m^v — 
938 MeV) the mass of the pion (nucleon) ^ . Higher-order corrections have been calculated for 
all dimension-six sources in Refs. |22t 125] . but they do not change the results significantly. The 
main difficulty (and uncertainty) in the calculation stems from the estimation of the low-energy 



constants (LECs) go and do — di in terms of the couplings in Eq. (61). Here we follow Ref. 



[m [25] and use naive dimensional analysis (NDA) ^26n27j. Hopefully, these estimates will be 
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Table 5: Bounds on the coupling constants of various dimension-six operators at Mqqd and Mw in 
units of (100TeV)-2. 





At = A/qcd a* = Mw 


du, d (/i) 
du, d (m) 
dw (m) 
ImSi (At) 
Im Es (At) 
ImELRi (a*) 
ImELRs M 


< {7.1, 6.8} < {15, 14} 

< {17, 8.0} < {18, 8.7} 

< 0.19 < 0.42 

< 2.3 < 0.51 

< 2.3 < 2.8 

< 2.3 < 1.7 

< 2.3 < 0.85 



replaced by lattice-QCD calculations. All estimates below are in terms of the dimension-six 
coupling constants at the scale Mqqd- 

For the quark EDM operators in Eq. (61) the main contribution to the neutron EDM comes 
from the short-range contribution do — di. The reason being that the non-electromagnetic 
pion-nucleon interaction is suppressed by the necessity of integrating out the photon appearing 
in the quark EDM operator. NDA gives do — di = 0{emuQudu, emdQddd) consistent with 
results obtained in quark models [H]. For the quark masses we use the values at = Mqcdj 
mu{MQCD) = 3.1 MeV and mrf(MQCD) = 6.5 MeV [38j. 

The quark CEDM operators induce both and, by including PT-even electromagnetic ef- 
fects, do — di. Both contributions in Eq. (63) appear formally at the same order (one expects, 
a priori, no cancellation among them) but the chiral logarithm, log(m^/m^) ~ —4, somewhat 
enhances the go contribution. NDA gives go = 0{mqdqMQQY)/ At:) consistent with a calculation 
using QCD sum rules [45j . 

The gluon CEDM conserves chiral symmetry and, as a consequence, its contribution to go 
(the LEG of a chiral-symmetry-breaking interaction) is suppressed by rUq/MqcD [Ml [25]. The 
main contribution to the nEDM stems from the short-range contribution which is estimated by 
do — di = 0{edwMQCD/4:Tr)- This estimate is approximately twice as large as a calculation based 
on QCD sum rules [46] which reflects the intrinsic uncertainty in these kind of estimations. 

Next are the FQPS operators with coupling constants ImSi^g. As pointed out in Refs. 
|241 [25] these operators conserve, just as the gluon CEDM, chiral symmetry such that for 
these sources is suppressed as well. The short-range terms are estimated as do — di = 
0(e(ImSi,8)MQCD/(47r)2). 

Finally, the FQLR operators do break chiral symmetry and, similar to the quark CEDM, 
the neutron EDM is determined by both terms in Eq. (63). However, for this dimension-six 
source is smaller than expected (something which goes beyond NDA) due to the smallness of 
the proton-neutron mass difference (for details, see Ref. [25]), such that the main contribution 
comes again from the short-range terms. NDA gives do — di = 0(e(ImHi_8)MQCD/(4vr)^). This 
estimate for the neutron EDM is about an order of magnitude smaller than the result obtained 
in Ref. [18]. A possible explanation for this discrepancy is given in Ref. [25]. To be on the safe 
side we will use the smaller estimate. 

The bounds on the dimension-six coupling constants in Eq. (61) are shown in Table [s] for 
the scales fi = Mqcd and fi = Mw- For the latter bounds, certain operators contribute to 
the neutron EDM in different ways. For example, a quark CEDM at Mw induces both a 
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Table 6: Bounds on the coupling constants of various dimension-six operators. All entries are dimen- 
sionless. 





= 1 TeV Mj- = 10 TeV 


{M^)dw (M;r) 
(M|.)ImEi {Mrf^) 
(M|.)ImE8 (M;r) 
(Af|.)ImHi (M;r) 
(Af|,)Imy'"''' (A'/;r) 

(M|,)6i' {Mrp) 


< {1.8, 1.8} ■ 10~^ < {2.1, 2.1} • 10"^ 

< {1.9, 0.91} • 10^=* < {1.7, 0.94} ■ 10"^ 

< 5.6 •10"'^ <7.0-10~^ 

< 3.2 -10"^ < 2.3 -10"^ 

< 3.3- 10"* < 2.4 -10"^ 

< 1.7 -10"* < 1.7 -10"^ 

< {8.9, 8.9} ■ 10-5 < {7.9, 7.9} • 10^^ 

< 2.4 -10"^ < 1.5- 10"^ 



quark EDM and CEDM at Mqcd which both contribute to the neutron EDM. In these cases we 
present the strongest bound and do not take into account possible cancellations between different 
contributions. The limits in Table [5] involve NDA factors and contain significant uncertainties. 
All bounds are of order (100 TeV)"^ apart from the bounds on the quark (C)EDMs which are 
approximately an order of magnitude weaker because they scale with the small quark mass. The 
running from the electroweak scale to low energies does not affect the bounds by a large amount. 
For most operators less than a factor two, the exception being the color-singlet FQPS operator 
whose bound is strengthened by a factor 5. The QCD corrections are thus not very significant, 
especially not in the light of the large uncertainties involved. 

It is also possible to set bounds on the couplings at energies above the electroweak scale. 
However, this depends on the particular value of in two ways. First, the coupling constants 
are expected to scale as ~ M^^. We therefore present bounds on the dimensionless combination 

Mj,Ci{M'^) where Ci denotes the various coupling constants. Second, the QCD corrections 
depend on the particular value of as well. In the first two columns of Table |6] we present 
bounds for two values of M^. The last three rows contain couplings which did not appear in Table 
[5] because the corresponding operators decoupled from the EFT below the electroweak scale. 
Although the gluon-Higgs operator only contributes to the qCEDM via a one-loop diagram, the 
bounds on 9' and dq are equally strong. 

It must be stressed that all bounds obtained here cannot simply be used to set bounds on 
M^. The reason is that the effective couplings can depend on dimensionless factors appearing 
in the high-energy theory and factors of 47r arising from integrating out the heavy fields. In 
fact, as argued in Ref. [17j, the couplings du,d du,d, dw, and 6' can only be generated at the loop 
level and are therefore associated with a factor l/(47r)^ ~ 10^^. It should be noted that these 
suppression factors are not general and do not always appear ^48j. For instance, compensation 
factors might arise if the fundamental theory is strongly coupled. In case such loop factors do 
appear, the bounds on the quark (C)EDMs are, even at the relatively low scale = 1 TeV, not 
particularly strong. The couplings Si,85 ^i, and Y'^''^ can be generated at tree level but they 
might still be suppressed by small couplings or phases. For example, in the minimal left-right 
symmetric model, Y'^'"^ is suppressed by the SM light-quark Yukawa couplings [ 39\ . 
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4 Suppressed operators 



In general, we expect the operators in Eq. (61) (which descend from Eq. (27)) to describe the 



dominant part of P and T violation due to physics beyond the SM in first-generation hadronic 
and nuclear systems. However, in specific models of new physics the associated coupling con- 
stants could be much smaller, or even zero, due to additional symmetry considerations or un- 
expected fine-tuning. In such models the most important operators might be operators which 
were neglected so far such as the dimension-six operators mentioned in Sec. [2] that obt ain ad- 
ditional suppression at low energies. Higher-dimensional operators might also become relevant, 
but we limit the discussion to the dimension-six operators. We assume that for some reason the 
operators of the last section are suppressed, thereby leaving room for the remaining dimension- 
six operators. In order to keep the discussion organized, we first discuss the quark-Higgs 



operators in Sec. 4.1 These operators give rise to a rich set of four-quark operators at low 



energies. In Sec. 4.2 we discuss the remaining operators which involve heavy gauge bosons. 



4.1 P- and T- violating quark-Higgs interactions 



We begin the analysis with Eq. (20) describing interactions between light quarks and Higgs 
bosons 

Cem^^v = hiv + h)(^v + ^^ (^i Im y'" + i Im Y"^ d-f^d) 

= ^ CgH OgH + . . . , (64) 

q=u,d 



where the dots denote terms with more than one Higgs boson and Eqs. (27) and (28) were 
applied. These interactions have already been considered in the previous section where they 
induced relatively large corrections to the quark CEDMs via two-loop diagrams involving top 
quarks [37]. Here we assume that these corrections, together with the other contributions to the 
quark (C)EDMs, are suppressed. This requires fine-tuning between different coupling constants 
which is not expected on general grounds but might occur in specific models of new physics. 

Starting at the high-energy scale M^, the Higgs-quark interactions are first evolved down to 
-^EW = ■ The anomalous dimension is given in Eq. ( |36| ) and the solution of the RGE is 

= { "4t^) { " C,„iMr)- (65) 



as{mt) J \as{Mf) J 

Numerically this becomes for two explicit values of 

CgH(MH/) = 1.2Cg^^(lTeV), CgH(MM/) = 1.3 (10 TeV). (66) 

Below Myi/ the Higgs boson is integrated out which at tree-level generates four-quark oper- 
ators among light quarks of the following form 

,qqqq — " ^72" ( '^"^"^ ("^W^) ^^^^^ (^^) ^')'^^ '^'^ 

+mdCuH{Mw) u-f^u dd + muCduiMw) uu dj^d] . (67) 
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Figure 4: One-loop diagrams contributing to the quark electric and chromo-electric dipole mo- 
ments. A square marks a interaction from Eq. (64). Other notation is as in Figs, [l] and [2} 
For simplicity only one possible ordering is shown here. 



Since Cqu = 0{l/M^) these interactions scale as 0{vmq/M'^M^) ~ 0{mq/vM^) such that 
they are not only suppressed by two powers of Mj' but also by the ratio of the light-quark mass 
to the Higgs vev, v ~ 246 GeV. The structure of the four-quark operators appearing in Eq. (67) 
is different from that in the previous section, which calls for an extension of the operator basis 
in Eqs. (27) and (47) with six additional four-quark operators 

■,2 



Ops, 
Ops, 



[uu d'jr^d + u^^u dd + du wy^d + dj^u ud), 

[ut^u dj^f^d + uj^t^-u df'd + dfu u-f^t^-d + dj^t^-u ut"-d) , 



O. 
O, 



4(ji 

'4'?8 



ids (qI^QQQ), 



(68) 



where again q G {u, d}. Using the Fierz identities 

- 1 ~ 1 

Ws ^7^"" dd = - {Ops, + Op5, ) + - OlRi + OlRs , 

Ws uudj^d = ^(0p5i + Op5i) - ^Olr, - OlRs, 
the coupling constants Ci{fi) at the scale are found to be 



(69) 



CpsAMw) 

Clr^Mw) 
C4,,(M-) 



CpsAM^y] 
1 



1 



1 



6 



Clpo(M, 



w) 



777 11 

C„h(M+) + -^Qh(M+) 



Ml 

ninV 



1 



CpsM 



<7f (M+) Ml 
Cps,{M-) = C,q,{M-) 



CgH(M+), 

0. (70) 



The quark-Higgs interactions induce the quark EDM and CEDM through one-loop dia- 
grams shown in Fig. |4j These diagrams are finite and give the following contributions 



Cg(M, 



w) 



-CM 



w) 



1 mqV 



3 ?Ti \ 



(71) 



which are smaller by a factor 0(10 ) compared to Eq. (40). This is approximately the level 
of fine-tuning needed in order to make the four-quark operators in Eq. (67) significant at low 
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energies. Eqs. (70) and (71) represent matching conditions at fi = Mw, meaning that the quark 
masses appearing there should be evaluated at this scale as well. The one-loop QCD running 
gives mq{M]y) ~ 0.58 mg(MQCD)- In all results below we will use niq to denote the quark mass 
at Mqcd- 

The next step is to evolve the operators down to Mqcd- We choose the basis 



CifJ-) = {Cq{fJ,),Cg{n),CwitJ'),Cps^itJ-),CpSsitJ-),CLRAfJ')^CLRs{fJ-), 
CPSI (^), CpSs (/i), C4gi (/i), Qqg (/i))^, 



with the extended RGE 

Up to first order in a^, 7 can be written as 



(72) 



(73) 



7 



( Tdipole 


7mix 





7mix 


'749, mix 





IPS 

















ILR 

















IPS 





V 











74g 



(74) 



The top-left entries of this matrix are given in Eqs. (32), (33), (34), and (51). The other 
entries denote the running and mixing of the four-quark operators in Eq. (68). As is already 
clear from the form of 7, these additional four-quark operators do not mix with the FQPS and 
FQLR operators. Furthermore, Op5,i,8 do not mix with 04g,i,8 at the one-loop level. These 
considerations justify the block-diagonal form of 7, apart from 7mix! 7mix, and 74g,mix which give 
rise to contributions to the quark (C)EDMs from the various four-quark interactions. Explicit 
calculation shows that 7P5 is identical to 74^ 



(4-3Af)(Af^-l) 



IPS = liq 



+ /3o 



I N-2 
N 



(N-l)'^{l+N){2+N) 

m 

2C2(iV) + 2^ + /3o 



(75) 



Apart from a sign Ops,i,s niixes into the quark (C)EDMs the same way as Op5,i,8 such that 
7mix = — 7mix- The four-quark operators containing a single flavor contribute to the dipoles via 



74g,mix = 4 



-C2{N) 
C2{N) - f 




(76) 



where it is implied that the operator containing only up (down) quarks induces only the up 
(down) quark (C)EDM. The results for the anomalous dimension matrices agree with those in 
Refs. [THl HH] after a basis transformation. 

Following the same procedure as in Sec. |3.2| we now find, apart from Table [3j additional 
contributions, shown in Table [7| to the quark (C)EDMs. The results for the gluon CEDM, the 
FQPS, and the FQLR operators are unchanged. The following results for the quark (C)EDMs 
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Table 7: Dependence of the quark electric and chromo-electric dipole moments and four-quark 
operators at Mqcd on the same four-quark operators at Mw. 





CpsAMw) 


Cps,{Mw) 






C,(Mqcd) 


f 0.070^ - 0.027) ^ 

V Qq J rnq 


(0.045^+0.0058) ^ 


0.085 


0.10 




-0.080^ 


0.018^ 

rriq 


-0.16 


0.037 


CpSiiMqcD) 


0.63 


-0.14 






Cpss{Mqcd) 


-0.063 


0.18 






C4,i(Mqcd) 






0.63 


-0.14 


C458(Mqcd) 






-0.063 


0.18 



are obtained for = 1 TeV 



dg(MQCD) 



dg(MQCD) 



rUn 



0.0020 ( 0.21 + In 

Mh 



0.00035 



0.15 
-0.028- 



0.020 ( 0.21 



0.0028 



0.038 

Ml 



Ml 



Imy'«'(lTeV) 



Imy"^(lTeV) + ^Imy'"(lTeV) 



rudV^ 



Ml 



Mh 



Ml 



Imy"?(lTeV) 



""^Imy'^ClTeV) + ^Imy'"(lTeV)j 



(77) 



H 



while for M^ = 10 TeV the results are trivially obtained from Eq. (66). Because the logarithm 
is rather large, Inniq/MH ~ —10, the one-loop diagrams contribute at the same order as the 
mixing of the four-quark operators into the dipoles. However, all terms receive a suppression of 
niq/v compared to Eqs. (54) and (55). 

The contributions to the various four-quark operators are suppressed by niq/v as well. For 
the four-quark operators appearing in Sec. |3} we find 



{ImSi(MQCD), ImSslMqco)} = {-1-6, 0.73} 
{ImHiR,(MQCD), ImHLR3(MQCD)} = -{0.34, 2.0} 



^Imy'"(lTeV) + ^lmy'^(lTeV) 
Ml ^ ' Ml ^ \ 



Ml 



Imy'"(lTeV) 



Ml 



Imy"^(lTeV) 



.(78) 



A comparison of Eqs. (77) and (78) makes it clear that the four-quark operators are larger by 
one to two orders of magnitude than the quark (C)EDMs. In the specific scenario discussed 
here, it seems safe to neglect the latter. 

Finally, there are the remaining six four-quark operators in Eq. ( 68 ) . In order to compare with 
Eq. (78) it is useful to redefine these operators without the factor and define new coupling 
constants 



^i,8(/^) = 5s(/^)C'pSi,8(m), 

%M = ffs(/^)C4g8(/^), 



(79) 
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for which we find 



{El (Mqcd) , Ss (Mqcd)} 
{Qg^i (Mqcd) , ^,,8 (AfqcD)} 



{-0.64, 0.064} 



Imy'"(l TeV) + -^lmY"^{l TeV) 



{-1.3, 0.13} -^Imr^(l TeV). 



(80) 



These results indicate that the color-singlet operators with couplings Si and ilgj , are larger by 
an order of magnitude than their color-octet counterparts Sg and ilgg. Thus, at low energies, it 
is sufficient to consider only the color-singlet operators which are of the same order as ImSi^g 
and IirElris- 

To summarize this section we write the Lagrangian at the scale Mqcd in case the dominant 
dimension-six operators at the electroweak scale are the Higgs-quark interactions in Eq. (64). 



Neglecting, as discussed above, the quark (C)EDMs and the Ss and Qq^s four-quark operators, 
the Lagrangian at Mqcd is given by 



6, set 2 



+i 



Im S 



2 

Im Eg 
2 



(uu dj^d + M75ii dd — du wy^d — d^^u ud) 

[ut^u dj^f^d + u-f^t^-u df'd - df'u u-f^t^-d - dj^t^-u uf^d) 



+ilmELRs {uR-f''t''dRdL-ff,t''uL - dRYt''uRULlt,t''dL) 

Si - - - - 

+i—{uu d'y^d + u^^u dd + du u'y^d + d'y^u ud) 

+i^}u^i{uj^uuu) + ind,i{d'l^ddd), 



(81) 



where the coupling constants are all evaluated at Mqcd and their values in term of the coupling 
constants in Eq. (64) are given in Eqs. (78) and (80). All coupling constants in Eq. (81) scale as 
{mq/v){l/Mj?j and, in general, are much smaller than the constants appearing in Eq. (61). We 



stress again that the new four-quark operators appearing in this section are only important in 
specific scenarios involving additional symmetries and/or fine-tuning. In the scenario sketched 
here there are seven relevant four-quark operators around Mqcd with only two independent 
coupling constants. 



4.2 P and T violation involving heavy gauge bosons 

In this section we deal with the remaining operators in Sec. [2] which contain heavy W^^-, 
and Higgs-boson fields. These operators will no longer exist in the effective theory below the 
electroweak scale, but they will contribute to operators containing light fields only. The operators 
in this section only give rise to low-energy PT violation via electroweak one-loop diagrams. The 
operators studied are the weak dipole moments of the quarks in Eqs. (10) and ( |12[ ) and the 
interactions among electroweak gauge and Higgs bosons in Eqs. (22), ([24), and (|25[). 

These operators contribute to quark (C)EDMs through various diagrams, however, in all cases 
the contributions are suppressed by electroweak coupling constants in the typical combination 
otw = e^/47r. The quark (C)EDMs are generated when the effective operators are evolved from 
M'p to Mew- Below Mew the heavy bosons can be integrated out and the running of the 
quark (C)EDMs to Mqcd can simply be obtained from the results in Sec. Is] In principle, the 
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Figure 5: One-loop diagrams contributing to the quark electric moment. The double wavy lines 
denote the propagation of and Z bosons. The black square denotes a vertex from one of 
the operators in Eq. (86). The other notation is as in Figs. [T] and [2j For simplicity only one 



possible ordering is shown here. 



quark (C)EDMs originate not only from the running from to Mewi but also from threshold 
corrections at Mw, i.e. the finite parts of the Feynman diagrams. Here we do not calculate 
these parts unless the operator in question does not produce a quark (C)EDM otherwise. This 
is only necessary for the interactions among gauge bosons with coupling constant in Eq. (22). 



For the other operators we assume the sizes of the induced quark (C)EDMs to be saturated by 
the running part. We do not expect this approximation to significantly alter the results. 

To calculate the running of the operators we will take the 0{aw) contributions as a pertur- 
bation to the QCD RGE. We only consider the electroweak corrections which induce a quark 
(C)EDM. That is, we do not consider the electroweak running of aw itself, nor that of the op- 
erators involving the heavy bosons. These effects would give rise to ©(a^) corrections to the 
induced quark (C)EDMs. The RGE is written as 



dC 
d In fj, 



7C, C = {Cg,C„CH? 



(82) 



where Ch stands for the coupling constant of one of the heavy-boson operators. These will be 
treated one at a time because their mixing can be neglected at the order we work. We write 
C = (7(°) + and 7 = 7(0) +7(1) where the superscript (1) denotes terms of ©(a^ ) {0{al,)). 
Up to 0{ai 



^^j, Eq. (82) splits into two equations 



9(7(0) 



dlnfi ^ ' 9 In ^ 

where 7^^^ describes the effects of QCD corrections, while 7*-^) describes the mixing of the 
operators due to electroweak corrections. The LO matrix can, for our purposes, be written as 



(0) 



^(0) c(i)+7(i) CW, 



(83) 



7 



(0) 



47r 



' 7e 7ge 
7, 
7H, 



(84) 



where Je,q,qe can be read off from Eq. ( |32[ ) (7e = — 7qe = 8C2(A^) and 7^ = 16C2(A'') — 4A'") 
while 7// and the matrix 7^^) depend on the heavy-boson operator under investigation. The 
electroweak effects can, at the order we work, be described by 



^(1) = ^ 
' 47r 



(0 















Vo 








(85) 



25 



What remains is the calculation of the unknown entries Tis^ cind 723^ for each operator. 
We begin by defining the relevant operators, 



Ob 
Ow 

OwB 
Odw 
Owu 

Owd 

Oza 



2 — F^^Za/s ) vh + 



FnuFap + 2 F^yZap ) vh + 

Sill 



FfiuZaf} 



h 

V 



--t^W;W~F^p) + ... 



9 



rud 



(uLicT^^'dii) W+ + i^idia^^'j^d) W^W- + h.c. 



9_ 
'%/2 

-^mqq^a^''-i^qZ^,^, 



i^{uia^-^^u) W^W- + h.c. 
9 



V2 



where the dots denote terms which are not relevant at 0{a^ 
scale Mj' are 



(86) 

The coupling constants at the 



CwB{Mf) -- 
Cwq{Mrp) -- 



(M^), 



Cdy^{Mrp) = -S^d^{Mrp), 



(87) 



The anomalous dimensions for the heavy boson operators are collected in Table [8j The 
electroweak loops for OwB-, Ody^, Ozq, and Owq were evaluated in Ref. |T2] in a different 
regularization scheme. However, for Owq not all diagrams were taken into account. 

The operators Ob, Ow, Owb, and induce a quark EDM through the diagrams in Fig. 
^ The first two operators contribute through Diagrams [5]l^a,b) and Ody^^ through Diagram [5|c). 
Owb contributes via all diagrams. None of these operators induces a quark CEDM at the 
one-loop level. As can be seen from Table [8| the anomalous dimensions of Ody^r vanish. The 



Table 8: The anomalous dimensions of the various heavy boson operators. Here stands for the third 
component of weak isospin of the external quark, i.e. 1/2 (—1/2) for the up (down) quark. 
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Owb 
Od„ 

Owq 
Ozq 



-4- 



-16 (1 - 
-16fl + 
I- 4s 



Tl/2~slC 

2 ' cos 20, 












g2 * ud 
,Tf-2Q^sl 










8C2(iV) 
8C2(iV) 
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(a) (b) (c) (d) (e) 



Figure 6: One-loop diagrams contributing to the quark electric and chromo-electric dipole mo- 
ments. The notation is as in Fig. [5] For simplicity only one possible ordering is shown here. 



Table 9: The contributions of the various heavy boson operators at to the q(C)EDMs at Mqcd in 
units of 10~^. A "— " indicates that there is no contribution at this order. 





Cb (ITeV) 


Cw (ITeV) 


Cws (ITeV) 


Cd„,(lTeV) 


CH'.ClTeV) 


Cz5(lTeV) 


Cu 


d (Mw) 


{1.8, 0.73} 


{3.6, 7.3} 


{6.2, 11} 


{0.65, 1.3} 


{-13, -15} 


{-2.4, 4.3} 


Cu 


d (Mw) 










{-2.2,-2.2} 


{-0.49, 0.88} 




Cb (lOTeV) 


Cw (lOTeV) 


C„.s( 10 TeV) 


Cd„,(10TeV) 


Cw',(10TeV) 


Cz,(10TeV) 


Cu 


d (Mw) 


{3.2, 1.3} 


{6.4, 13} 


{11, 19} 


{0.65, 1.3} 


{-25, -29} 


{-4.9, 8.8} 


Cu 


d (Mw) 










{-3.6,-3.6} 


{-0.81, 1.5} 



reason is that for this operator Diagram [sjj^c) is finite in dimensional regularization. However, 
the diagram does have a finite part which was evaluated in several regularization schemes in Ref. 
|49j . where it was found to be scheme dependent. We will use the result found in dimensional 
regularization. In our notation the operator Od^ induces the quark EDM |49] 



C„,.(Mi.) = ±^^Q„„ (88) 

where the plus (minus) sign is for the up (down) quark EDM. 

The weak dipole moments, Ozq and Owq-, contribute to four-quark operators at tree level. 
However, this necessarily involves one power of the exchanged momentum because of the deriva- 
tive acting on the gauge fields. The coupling constants of these effective dimension-seven op- 
erators then scale as m^/M^ and are heavily suppressed. The same holds for one-loop 
contributions to four-quark operators. Larger effects come from one-loop diagrams shown in Fig. 
[6]j^a-d) and [oj^e) contributing to, respectively, light-quark EDMs and CEDMs. The weak dipole 
moments are the only operators in this section which produce a quark CEDM. Furthermore, 
they are the only operators which are affected by QCD corrections and, hence, have a nonzero 

Employing all the anomalous dimensions, the finite contribution from O^^, and using |38) 

a^{Mw) ^ sl{Mw) ^ 0.23, (89) 

we calculate the induced quark (C)EDMs at M\\r. The results for = 1, lOTeV are given in 
Table [9j The quark EDM entries for Cb-, Cw-, and CwB for = lOTeV are almost twice as 
large as those for = 1 TeV. This can be understood from the logarithmic dependence on 
Mrp, while, simultaneously, QCD corrections suppress the contributions. The same holds for the 
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Table 10: Bounds on the couplings of various operators involving heavy bosons. All entries are 
dimensionless. 





Mf^l TeV Mf = 10 TeV 


{M1)Cb (M;r) 
{mI)Cw {Mj.) 
{m{)Cwb {Mf) 

{Mj,)Cwu,d{M^) 


< 8.1 -10"^ <4.6 

< 1.9 10"^ < 1.1 

< 1.3 -10"^ <0.74 
<0.11 < 11 

< {1.0, 0.84} ■ 10-2 < {0.53, 0.45} 

< {5.3, 2.8} • 10^2 < {2.7, 1.4} 



quark CEDM entries for Cwq and Czq ■ The quark EDM entries for these last two operators are 
more complicated because they get additional contributions when the induced quark CEDM is 
evolved to lower energies. 

The quark (C)EDMs can be run down to Mqcd using results from Section |3j Following 
Section 3.3 we use the neutron EDM limit to set bounds on the couplings of the heavy-boson 



operators. Using the same estimates as before, we obtain the results in Table [T0| 

The limits on the different operators are all of similar size apart from those on which are 
somewhat weaker. A comparison with Table [6] shows that the bounds on the heavy-boson oper- 
ators are approximately an order of magnitude weaker than the bounds on the quark (C)EDMs. 
While they are a factor of 0(10^-10^) weaker than the limits on Iml"''', ImHi, and ImSi^g. We 
mention again that these bounds cannot simply be used to constrain M^. The effective couplings 
will depend on dimensionless quantities coming from the high-energy theory. In fact, in theories 
for which the results of Ref. [U] hold all operators discussed in this section cannot be gener- 
ated at tree level and so their couplings will be suppressed by loop factors ~ l/(47r)^ ~ lO"^. 
Since the bounds at 1 TeV, let alone 10 TeV, are all of that level or weaker, for such theories 
the neutron EDM does not significantly constrain the coupling constants of the dimension-six 
operators at this scale. 

Some of the operators studied here have been under recent investigation as they can modify 
the /i — )■ 77 rate [501 [361 H] . lu this context, Ref. [51] used the electron EDM limit to set bounds 
on Ow-, Ob, Owe, and Od^y- (Note that the obtained limits are, for all four operators, stronger 
than limits obtained from accelerator-based experiments [52]). The obtained constraints are 
approximately an order of magnitude stronger than those derived here. The main difference is 
due to the fact that the electron EDM limit ^53j is about 30 times stronger than the neutron 
EDM limit, although this effect is softened by the smallness of the electron mass with respect 
to the light-quark mass. Furthermore, the electron EDM is not subject to QCD corrections 
which suppress the quark EDMs by a small amount when they are evolved to lower energies. 
Overall, the electron and neutron EDM searches are complementary because they are sensitive 
to different sources. For example, if the dominant source is the operator Ow the electron 
and neutron EDM would be of approximately the same size, whereas if, say, the gluon CEDM 
is the dominant source, the neutron EDM is expected to be much larger. This illustrates that 
measurements on different systems are needed to disentangle the fundamental mechanism. 
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5 Discussion and conclusion 



A measurement of a nonzero EDM in any of the upcoming experiments would be a major 
breakthrough. Within the SM, a non-zero hadronic or nuclear EDM at the current experimental 
accuracies, can only originate in the 9 term because P and T violation from the quark- mixing 
matrix is simply too small. However, the severe suppression of 9 leaves room for physics beyond 
the SM. The effects of new physics can be parametrized by effective higher-dimensional operators 
which start at dimension six [10]. In this article we have performed a systematic study of the 
dimension-six operators relevant to flavor-diagonal P and T violation in hadronic and nuclear 
systems. In particular, we have investigated the evolution of these operators from the energy 
scale where they originate (M^), assumed to be significantly larger than the electroweak scale 
(Mew), to the QCD scale (Mqcd ~ IGeV). 

The effective operators appearing at contain SM fields only and obey the SM SUdS) x 
SUl{2) X f/y(l) gauge symmetries. The couplings of these operators are all proportional to two 
inverse powers of Mj'. Gauge invariance puts strong constraints on the form of the effective 
operators. For example, only two four-quark operators containing first-generation quarks, 
i.e. the FQPS operators, are allowed. Operators such as the quark EDM and CEDM which 
have canonical dimension five at low energies are actually dimension-six operators in disguise. 
Their chiral-symmetry-breaking properties forces them to be coupled to the Higgs field at high 
energies in order to preserve the SM gauge symmetries. Although at lower energies the Higgs 
field takes on its vev, the operators still scale as [12] • 

In order to study the effects of the dimension-six operators on low-energy observables, such 
as EDMs, it is necessary to evolve the operators to these lower-energy scales. In the process 
heavy SM fields decouple and must be integrated out. Simultaneously, it is necessary to include 
the effects of QCD and, in some cases, electroweak renormalization-group running. While doing 
so, we have found it convenient to divide the dimension-six operators at into two sets. The 
division is based on the size of the induced low-energy operators. The first set contains the 
dimension-six operators which induce low-energy operators that also scale as M^^. That is, 
the operators obtain no additional suppression when evolved to lower energies apart from minor 
suppressions due to QCD corrections. The second set contains the remaining operators which 
are suppressed by SM factors when run to Mqcd- Of course, EDM limits put the strongest 
constraints on operators in the first set. 



Set 1: Operators in the first set have been discussed in Sec. [3} The set consists of operators 
which contain only light fields such as the gluon CEDM or the FQPS operators. The shape of 
these operators stays the same all the way from M^ to Mqcd and, apart from QCD corrections, 
their scalings are not altered. Other operators in the first set are the quark EDMs and CEDMs. 
At high energies these operators couple to a Higgs field, but these can be replaced by their vev 
and there are no suppression factors at low energies. The next operator describes a derivative 
interaction among right-handed quarks and two Higgs fields (see Eq. (17)). When the Higgs field 
takes on its vev, an interaction among right-handed quarks and a W boson remains. Integrating 
out the W boson induces an additional four-quark operator (the FQLR operator) [33] . The 
W exchange suppresses the operator by M^, but this is compensated by the two powers of 
V appearing in the operator. Finally, the first set contains quark-Higgs and gluon-Higgs 
interactions. When the Higgs field takes its vev, the resulting terms renormalize, respectively, 
the SM Yukawa couplings and the 9 term. Via loop diagrams the quark-Higgs [37] and gluon- 
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Higgs interactions induces contributions to the quark (C)EDM. The associated loop suppressions 
are not very stringent and we therefore keep these operators in the first set. 



The low-energy form of the operators in the first set has been derived in Sec. 3.2 Around 



Mqcd there are seven relevant operators among light quarks, gluons, and photons. They 
consist of the quark EDM, the quark and gluon CEDM, and four four-quark operators. The 
first three of these operators have been extensively considered in the literature. The four-quark 
operators are less often taken into account, possibly due to the belief that they should be 
proportional to the quark masses and are therefore effectively higher-dimensional operators in 
disguise. Although this arguments holds for some four-quark operators it does not hold for the 
FQPS and FQLR operators which genuinely scale as M^'^ and are in general not suppressed 
compared to the dipole operators. Model-independent studies of hadronic and nuclear EDMs 
due to beyond-the-SM physics should use the seven operators at Mqcd as a starting point. Only 
model-dependent statements can single out particular operators in this set. 

The FQPS and FQLR operators have very different origins. As mentioned, the FQPS oper- 
ators conserve the SM gauge symmetries and appear directly at M^. There is, a 'priori^ no link 
between the coupling constants of the color-singlet and color-octet FQPS operators since both 
are independently allowed by the gauge symmetries. Of course, such a link can exist in specific 
models of new physics. QCD corrections particularly affect the color-singlet operator when it is 
evolved to Mqcd as can be seen from Eqs. (57) and (58). The FQLR operators do not conserve 
the SU2{L) gauge symmetry and are therefore not allowed at M^. Nevertheless, the color-singlet 
FQLR is generated below the electroweak scale after a VK-boson exchange [S^. The color-octet 
FQLR operator is only induced when its color-singlet partner is evolved to lower energies. The 
two operators are therefore not independent and both depend on the same coupling constant 
Si appearing in Eq. (17) [25]. Eq. (60) shows that at Mqcd the coefficient of the color-octet 
FQLR operator is even slightly larger than that of the color-singlet operator. 

In Sec. 3.3 we have used the neutron EDM limit to set bounds on the seven operators at 



Mqcd using xPT calculations [HI [25]. By use of the results in Sees. 3.1 and 3.2 these limits 
are translated to limits on all dimension-six operators in the first set at energies around M^. 
This makes it possible to quickly obtain restrictions imposed by the neutron EDM on specific 
high-energy models. After performing a matching calculation between the high-energy model 
and the effective dimension-six operators, these restrictions can be read off immediately. 

There are interesting proposals to measure the EDMs of light nuclei in storage rings with high 
accuracy [9j. Such an experimental program could give important complementary information 
on the dominant mechanism. In Refs. j2H a strategy has been proposed to 

(partially) disentangle the/^ sources at the energy scale Mqcd- For example, a deuteron EDM 
significantly larger than the sum of the neutron and the proton EDM would point towards a 
quark CEDM or the FQLR operators ^20t i24j . On the other hand, if the quark EDM is dominant 
one would expect the deuteron EDM to be close to the sum of the nucleon EDMs [21]. This 
strategy combined with the results obtained here can be used to (partially) disentangle the 
various sources at high energies, if any exist. 



Set 2: The second set of dimension-six operators at M^ consists of the remaining operators 
which do suffer from additional suppression when evolved to lower energies. We first discuss 
quark-Higgs interactions which already appeared in the first set because they induced relatively 
large corrections to the quark CEDM via loops. At tree- level the quark-Higgs couplings induce 
four-quark operators which are different from the FQPS and FQLR operators. However, 
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the couplings of these operators scale as {mq/v)M^ and are thus additionally suppressed by 
the ratio of the light-quark mass to the Higgs vev. It is interesting that when the quark-Higgs 
operators are evolved to energies around Mqcd, all ten possible/^ four-quark operators among 
up and down quarks are induced, albeit with only two independent couplings and the coefficients 
of three operators are an order of magnitude smaller than the others. The QCD evolution of 
these operators was studied in a different basis in Refs. [ TSj I19j . The set contains the FPQS 
and FQLR operators but also six new structures. However, because the quark-Higgs operators 
induce a larger contribution to the quark CEDM, the additional six four-quark operators can in 
general be neglected. 

The other dimension-six operators in the second set describe interactions involving elec- 
troweak gauge and/or Higgs bosons. The operators consist of the quark weak dipole moments 
and interactions among three gauge bosons and/or two gauge bosons and a Higgs boson. 
The quark weak dipole moments are similar to the quark EDM but the photon is replaced by 
a W or Z boson. In fact, as seen from Eqs. (11), and (13) these operators are related by 
gauge symmetry. The purely bosonic interactions originate in four independent gauge-invariant 
operators flO\. 

All operators contain at least one heavy field which decouples below the electroweak scale 
generating contributions to the quark EDMs and CEDMs via electroweak one-loop diagrams 
[I2] . These contributions are typically suppressed by au;/47r. At the one-loop level there are 
no dimension-six four-quark operators generated. In Sec. |4.2| we have calculated the induced 
quark (C)EDMs and used the neutron EDM limit to set bounds on the various dimension-six 
operators. These bounds are weaker than the bounds on the operators in the first set, but by 
how much depends on the particular operator. 

The purely bosonic dimension-six operators induce, apart from quark EDMs, also contribu- 
tions to the electron EDM [51j via very similar diagrams. The limits obtained from the electron 
EDM are stronger than from the neutron EDM (both are stronger than limits from accelerator 
experiments |5H I52j) but this might change depending on which EDM limit is improved first. 
What is interesting is that the bosonic operators induce an electron and neutron EDM of similar 
size, while the other dimension-six operators discussed in this paper induce a significantly larger 
neutron EDM. Observation of such a pattern in future experiments could provide a hint towards 
the fundamental source of P and T violation. 



Finally, we should say that our set of dimension-six operators around is far from com- 
plete since we did not include operators involving heavier quarks. Allowing such fields would 
increase the number of operators significantly, especially if one allows for generation-changing 
operators. Most importantly, we did not consider operators involving strange quarks which are 
still present in the EFT at Mqcd and can have important consequences for hadronic and nu- 
clear EDMs. We leave a study of operators involving strangeness to future work. In general, 
we expect operators involving heavier quarks than up, down, and strange to give suppressed 
contributions to operators around Mqcd due to the need of integrating out the heavier quarks. 
However, exceptions to this argument exist, see, for example, discussions in Refs. [Tj I19j. 

In conclusion, we have performed a systematic study of parity- and time-reversal violating 
operators of dimension six which originate in physics beyond the SM. Beginning at the high- 
energy scale where these operators originate and their forms are constrained by gauge invariance, 
we have evolved the operators down to the electroweak scale and subsequently to hadronic scales. 
We have derived a set of operators which is expected to dominate hadronic and nuclear EDMs 



31 



due to physics beyond the SM. Furthermore, we have obtained quantitative relations between 
these operators and the original dimension-six operators at the high-energy scale. 
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